Quantitative Evaluation of Bone Microstructure using High-Resolution Extremity Cone-Beam CT with a CMOS Detector.

PURPOSE: A high-resolution cone-beam CT (CBCT) system for extremity imaging has been developed using a custom complementary metal-oxide-semiconductor (CMOS) x-ray detector. The system has spatial resolution capability beyond that of recently introduced clinical orthopedic CBCT. We evaluate performance of this new scanner in quantifying trabecular microstructure in subchondral bone of the knee. METHODS: The high-resolution scanner uses the same mechanical platform as the commercially available Carestream OnSight 3D extremity CBCT, but replaces the conventional amorphous silicon flat-panel detector (a-Si:H FPD with 0.137 mm pixels and a ~0.7 mm thick scintillator) with a Dalsa Xineos3030 CMOS detector (0.1 mm pixels and a custom 0.4 mm scintillator). The CMOS system demonstrates ~40% improved spatial resolution (FWHM of a ~0.1 mm tungsten wire) and ~4× faster scan time than FPD-based extremity CBCT (FPD-CBCT). To investigate potential benefits of this enhanced spatial resolution in quantitative assessment of bone microstructure, 26 trabecular core samples were obtained from four cadaveric tibias and imaged using FPD-CBCT (75 µm voxels), CMOS-CBCT (75 µm voxels), and reference micro-CT (µCT, 15 µm voxels). CBCT bone segmentations were obtained using local Bernsen's thresholding combined with global histogram-based pre-thresholding; µCT segmentation involved Otsu's method. Measurements of trabecular thickness (Tb.Th), spacing (Tb.Sp), number (Tb.N) and bone volume (BV/TV) were performed in registered regions of interest in the segmented CBCT and µCT reconstructions. RESULTS: CMOS-CBCT achieved noticeably improved delineation of trabecular detail compared to FPD-CBCT. Correlations with reference µCT for metrics of bone microstructure were better for CMOS-CBCT than FPD-CBCT, in particular for Tb.Th (increase in Pearson correlation from 0.84 with FPD-CBCT to 0.96 with CMOS-CBCT) and Tb.Sp (increase from 0.80 to 0.85). This improved quantitative performance of CMOS-CBCT is accompanied by a reduction in scan time, from ~60 sec for a clinical high resolution protocol on FPD-CBCT to ~17 sec for CMOS-CBCT. CONCLUSION: The CMOS-based extremity CBCT prototype achieves improved performance in quantification of bone microstructure, while retaining other diagnostic capabilities of its FPD-based precursor, including weight-bearing imaging. The new system offers a promising platform for quantitative imaging of skeletal health in osteoporosis and osteoarthritis.

Quantitative Cone-Beam CT of Bone Mineral Density Using Model-Based Reconstruction.

PURPOSE: We develop and validate a model-based framework for artifact correction and image reconstruction to enable application of Cone-Beam CT (CBCT) in quantitative assessment of bone mineral density (BMD). Compared to conventional quantitative CT, this approach does not require a BMD calibration phantom in the field-of-view during an object scan. METHODS: The quantitative CBCT (qCBCT) imaging framework combined fast Monte Carlo (MC) scatter estimation, accurate models of detector response, and polyenergetic Poisson likelihood (PolyPL, Elbakri et al 2003). The underlying object model assumed that the tissues were ideal mixtures of water and calcium carbonate (CaCO3). Accuracy and reproducibility of qCBCT was evaluated in benchtop test-retest studies emulating a compact extremity CBCT system (axis-detector distance=56 cm, 90 kVp x-ray beam, ~16 mGy central dose). Various arrangements of Ca inserts (50-500 mg/mL) were placed in water cylinders of ~11 cm to ~15 cm diameter and scanned at multiple positions inside the field-of-view for a total of 20 configurations. In addition, a cadaveric ankle was imaged in five configurations (with and without Ca inserts and water bath). Coefficient of variation (CV) of BMD values across different experimental configurations was used to assess reproducibility under varying imaging conditions. The performance of the model-based qCBCT framework (MC + PolyPL) was compared to FDK with water beam hardening correction and MC scatter correction. RESULTS: The PolyPL framework achieved accuracy of 20 mg/mL or better across all insert densities and experimental configurations. By comparison, the accuracy of the FDK-based BMD estimates deteriorated with higher mineralization, resulting in ~120 mg/mL error for a 500 mg/mL Ca insert. Additionally, the model-based approach mitigated residual streaks that were present in FDK reconstructions. The CV of both methods was ~15% at 50 mg/mL Ca and less than ~8% for higher density inserts, where the PolyPL framework achieved 20-25% lower CV than the FDK-based approach. CONCLUSION: Accurate and reproducible BMD measurements can be achieved in extremity CBCT, supporting clinical applications in quantitative monitoring of fracture risk, osteoporosis treatment, and early osteoarthritis.

Robust Quantitative Assessment of Trabecular Microarchitecture in Extremity Cone-Beam CT Using Optimized Segmentation Algorithms.

PURPOSE: In-vivo evaluation of bone microarchitecture remains challenging because of limited resolution of conventional orthopaedic imaging modalities. We investigate the performance of flat-panel detector extremity Cone-Beam CT (CBCT) in quantitative analysis of trabecular bone. To enable accurate morphometry of fine trabecular bone architecture, advanced CBCT pre-processing and segmentation algorithms are developed. METHODS: The study involved 35 transilliac bone biopsy samples imaged on extremity CBCT (voxel size 75 µm, imaging dose ~13 mGy) and gold standard µCT (voxel size 7.67 µm). CBCT image segmentation was performed using (i) global Otsu's thresholding, (ii) Bernsen's local thresholding, (iii) Bernsen's local thresholding with additional histogram-based global pre-thresholding, and (iv) the same as (iii) but combined with contrast enhancement using a Laplacian Pyramid. Correlations between extremity CBCT with the different segmentation algorithms and gold standard µCT were investigated for measurements of Bone Volume over Total Volume (BV/TV), Trabecular Thickness (Tb.Th), Trabecular Spacing (Tb.Sp), and Trabecular Number (Tb.N). RESULTS: The combination of local thresholding with global pre-thresholding and Laplacian contrast enhancement outperformed other CBCT segmentation methods. Using this optimal segmentation scheme, strong correlation between extremity CBCT and µCT was achieved, with Pearson coefficients of 0.93 for BV/TV, 0.89 for Tb.Th, 0.91 for Tb.Sp, and 0.88 for Tb.N (all results statistically significant). Compared to a simple global CBCT segmentation using Otsu's algorithm, the advanced segmentation method achieved ~20% improvement in the correlation coefficient for Tb.Th and ~50% improvement for Tb.Sp. CONCLUSIONS: Extremity CBCT combined with advanced image pre-processing and segmentation achieves high correlation with gold standard µCT in measurements of trabecular microstructure. This motivates ongoing development of clinical applications of extremity CBCT in in-vivo evaluation of bone health e.g. in early osteoarthritis and osteoporosis.

High-Resolution Extremity Cone-Beam CT with a CMOS Detector: Evaluation of a Clinical Prototype in Quantitative Assessment of Bone Microarchitecture.

PURPOSE: A prototype high-resolution extremity cone-beam CT (CBCT) system based on a CMOS detector was developed to support quantitative in vivo assessment of bone microarchitecture. We compare the performance of CMOS CBCT to an amorphous silicon (a-Si:H) FPD extremity CBCT in imaging of trabecular bone. METHODS: The prototype CMOS-based CBCT involves a DALSA Xineos3030 detector (99 µm pixels) with 400 µm-thick CsI scintillator and a compact 0.3 FS rotating anode x-ray source. We compare the performance of CMOS CBCT to an a-Si:H FPD scanner built on a similar gantry, but using a Varian PaxScan2530 detector with 0.137 mm pixels and a 0.5 FS stationary anode x-ray source. Experimental studies include measurements of Modulation Transfer Function (MTF) for the detectors and in 3D image reconstructions. Image quality in clinical scenarios is evaluated in scans of a cadaver ankle. Metrics of trabecular microarchitecture (BV/TV, Bone Volume/Total Volume, TbSp, Trabecular Spacing, and TbTh, trabecular thickness) are obtained in a human ulna using CMOS CBCT and a-Si:H FPD CBCT and compared to gold standard µCT. RESULTS: The CMOS detector achieves ~40% increase in the f20 value (frequency at which MTF reduces to 0.20) compared to the a-Si:H FPD. In the reconstruction domain, the FWHM of a 127 µm tungsten wire is also improved by ~40%. Reconstructions of a cadaveric ankle reveal enhanced modulation of trabecular structures with the CMOS detector and soft-tissue visibility that is similar to that of the a-Si:H FPD system. Correlations of the metrics of bone microarchitecture with gold-standard µCT are improved with CMOS CBCT: from 0.93 to 0.98 for BV/TV, from 0.49 to 0.74 for TbTh, and from 0.9 to 0.96 for TbSp. CONCLUSION: Adoption of a CMOS detector in extremity CBCT improved spatial resolution and enhanced performance in metrics of bone microarchitecture compared to a conventional a-Si:H FPD. The results support development of clinical applications of CMOS CBCT in quantitative imaging of bone health.

High-resolution extremity cone-beam CT with a CMOS detector: Task-based optimization of scintillator thickness.

PURPOSE: CMOS x-ray detectors offer small pixel sizes and low electronic noise that may support the development of novel high-resolution imaging applications of cone-beam CT (CBCT). We investigate the effects of CsI scintillator thickness on the performance of CMOS detectors in high resolution imaging tasks, in particular in quantitative imaging of bone microstructure in extremity CBCT. METHODS: A scintillator thickness-dependent cascaded systems model of CMOS x-ray detectors was developed. Detectability in low-, high- and ultra-high resolution imaging tasks (Gaussian with FWHM of ~250 µm, ~80 µm and ~40 µm, respectively) was studied as a function of scintillator thickness using the theoretical model. Experimental studies were performed on a CBCT test bench equipped with DALSA Xineos3030 CMOS detectors (99 µm pixels) with CsI scintillator thicknesses of 400 µm and 700 µm, and a 0.3 FS compact rotating anode x-ray source. The evaluation involved a radiographic resolution gauge (0.6-5.0 lp/mm), a 127 µm tungsten wire for assessment of 3D resolution, a contrast phantom with tissue-mimicking inserts, and an excised fragment of human tibia for visual assessment of fine trabecular detail. RESULTS: Experimental studies show ~35% improvement in the frequency of 50% MTF modulation when using the 400 µm scintillator compared to the standard nominal CsI thickness of 700 µm. Even though the high-frequency DQE of the two detectors is comparable, theoretical studies show a 14% to 28% increase in detectability index (d'2) of high- and ultrahigh resolution tasks, respectively, for the detector with 400 µm CsI compared to 700 µm CsI. Experiments confirm the theoretical findings, showing improvements with the adoption of 400 µm panel in the visibility of the radiographic pattern (2× improvement in peak-to-through distance at 4.6 lp/mm) and a 12.5% decrease in the FWHM of the tungsten wire. Reconstructions of the tibial plateau reveal enhanced visibility of trabecular structures with the CMOS detector with 400 µm scinitllator. CONCLUSION: Applications on CMOS detectors in high resolution CBCT imaging of trabecular bone will benefit from using a thinner scintillator than the current standard in general radiography. The results support the translation of the CMOS sensor with 400 µm CsI onto the clinical prototype of CMOS-based extremity CBCT.

Motion compensation in extremity cone-beam CT using a penalized image sharpness criterion.

Cone-beam CT (CBCT) for musculoskeletal imaging would benefit from a method to reduce the effects of involuntary patient motion. In particular, the continuing improvement in spatial resolution of CBCT may enable tasks such as quantitative assessment of bone microarchitecture (0.1 mm-0.2 mm detail size), where even subtle, sub-mm motion blur might be detrimental. We propose a purely image based motion compensation method that requires no fiducials, tracking hardware or prior images. A statistical optimization algorithm (CMA-ES) is used to estimate a motion trajectory that optimizes an objective function consisting of an image sharpness criterion augmented by a regularization term that encourages smooth motion trajectories. The objective function is evaluated using a volume of interest (VOI, e.g. a single bone and surrounding area) where the motion can be assumed to be rigid. More complex motions can be addressed by using multiple VOIs. Gradient variance was found to be a suitable sharpness metric for this application. The performance of the compensation algorithm was evaluated in simulated and experimental CBCT data, and in a clinical dataset. Motion-induced artifacts and blurring were significantly reduced across a broad range of motion amplitudes, from 0.5 mm to 10 mm. Structure similarity index (SSIM) against a static volume was used in the simulation studies to quantify the performance of the motion compensation. In studies with translational motion, the SSIM improved from 0.86 before compensation to 0.97 after compensation for 0.5 mm motion, from 0.8 to 0.94 for 2 mm motion and from 0.52 to 0.87 for 10 mm motion (~70% increase). Similar reduction of artifacts was observed in a benchtop experiment with controlled translational motion of an anthropomorphic hand phantom, where SSIM (against a reconstruction of a static phantom) improved from 0.3 to 0.8 for 10 mm motion. Application to a clinical dataset of a lower extremity showed dramatic reduction of streaks and improvement in delineation of tissue boundaries and trabecular structures throughout the whole volume. The proposed method will support new applications of extremity CBCT in areas where patient motion may not be sufficiently managed by immobilization, such as imaging under load and quantitative assessment of subchondral bone architecture.

Image-Based Motion Compensation for High-Resolution Extremities Cone-Beam CT.

PURPOSE: Cone-beam CT (CBCT) of the extremities provides high spatial resolution, but its quantitative accuracy may be challenged by involuntary sub-mm patient motion that cannot be eliminated with simple means of external immobilization. We investigate a two-step iterative motion compensation based on a multi-component metric of image sharpness. METHODS: Motion is considered with respect to locally rigid motion within a particular region of interest, and the method supports application to multiple locally rigid regions. Motion is estimated by maximizing a cost function with three components: a gradient metric encouraging image sharpness, an entropy term that favors high contrast and penalizes streaks, and a penalty term encouraging smooth motion. Motion compensation involved initial coarse estimation of gross motion followed by estimation of fine-scale displacements using high resolution reconstructions. The method was evaluated in simulations with synthetic motion (1-4 mm) applied to a wrist volume obtained on a CMOS-based CBCT testbench. Structural similarity index (SSIM) quantified the agreement between motion-compensated and static data. The algorithm was also tested on a motion contaminated patient scan from dedicated extremities CBCT. RESULTS: Excellent correction was achieved for the investigated range of displacements, indicated by good visual agreement with the static data. 10-15% improvement in SSIM was attained for 2-4 mm motions. The compensation was robust against increasing motion (4% decrease in SSIM across the investigated range, compared to 14% with no compensation). Consistent performance was achieved across a range of noise levels. Significant mitigation of artifacts was shown in patient data. CONCLUSION: The results indicate feasibility of image-based motion correction in extremities CBCT without the need for a priori motion models, external trackers, or fiducials.

Cone-Beam CT of Traumatic Brain Injury Using Statistical Reconstruction with a Post-Artifact-Correction Noise Model.

Traumatic brain injury (TBI) is a major cause of death and disability. The current front-line imaging modality for TBI detection is CT, which reliably detects intracranial hemorrhage (fresh blood contrast 30-50 HU, size down to 1 mm) in non-contrast-enhanced exams. Compared to CT, flat-panel detector (FPD) cone-beam CT (CBCT) systems offer lower cost, greater portability, and smaller footprint suitable for point-of-care deployment. We are developing FPD-CBCT to facilitate TBI detection at the point-of-care such as in emergent, ambulance, sports, and military applications. However, current FPD-CBCT systems generally face challenges in low-contrast, soft-tissue imaging. Model-based reconstruction can improve image quality in soft-tissue imaging compared to conventional filtered backprojection (FBP) by leveraging high-fidelity forward model and sophisticated regularization. In FPD-CBCT TBI imaging, measurement noise characteristics undergo substantial change following artifact correction, resulting in non-negligible noise amplification. In this work, we extend the penalized weighted least-squares (PWLS) image reconstruction to include the two dominant artifact corrections (scatter and beam hardening) in FPD-CBCT TBI imaging by correctly modeling the variance change following each correction. Experiments were performed on a CBCT test-bench using an anthropomorphic phantom emulating intra-parenchymal hemorrhage in acute TBI, and the proposed method demonstrated an improvement in blood-brain contrast-to-noise ratio (CNR = 14.2) compared to FBP (CNR = 9.6) and PWLS using conventional weights (CNR = 11.6) at fixed spatial resolution (1 mm edge-spread width at the target contrast). The results support the hypothesis that FPD-CBCT can fulfill the image quality requirements for reliable TBI detection, using high-fidelity artifact correction and statistical reconstruction with accurate post-artifact-correction noise models.

Dual-Energy Imaging of Bone Marrow Edema on a Dedicated Multi-Source Cone-Beam CT System for the Extremities.

PURPOSE: Arthritis and bone trauma are often accompanied by bone marrow edema (BME). BME is challenging to detect in CT due to the overlaying trabecular structure but can be visualized using dual-energy (DE) techniques to discriminate water and fat. We investigate the feasibility of DE imaging of BME on a dedicated flat-panel detector (FPD) extremities cone-beam CT (CBCT) with a unique x-ray tube with three longitudinally mounted sources. METHODS: Simulations involved a digital BME knee phantom imaged with a 60 kVp low-energy beam (LE) and 105 kVp high-energy beam (HE) (+0.25 mm Ag filter). Experiments were also performed on a test-bench with a Varian 4030CB FPD using the same beam energies as the simulation study. A three-source configuration was implemented with x-ray sources distributed along the longitudinal axis and DE CBCT acquisition in which the superior and inferior sources operate at HE (and collect half of the projection angles each) and the central source operates at LE. Three-source DE CBCT was compared to a double-scan, single-source orbit. Experiments were performed with a wrist phantom containing a 50 mg/ml densitometry insert submerged in alcohol (simulating fat) with drilled trabeculae down to ~1 mm to emulate the trabecular matrix. Reconstruction-based three-material decomposition of fat, soft tissue, and bone was performed. RESULTS: For a low-dose scan (36 mAs in the HE and LE data), DE CBCT achieved combined accuracy of ~0.80 for a pattern of BME spherical lesions ranging 2.5 - 10 mm diameter in the knee phantom. The accuracy increased to ~0.90 for a 360 mAs scan. Excellent DE discrimination of the base materials was achieved in the experiments. Approximately 80% of the alcohol (fat) voxels in the trabecular phantom was properly identified both for single and 3-source acquisitions, indicating the ability to detect edemous tissue (water-equivalent plastic in the body of the densitometry insert) from the fat inside the trabecular matrix (emulating normal trabecular bone with significant fraction of yellow marrow). CONCLUSION: Detection of BME and quantification of water and fat content were achieved in extremities DE CBCT with a longitudinal configuration of sources providing DE imaging in a single gantry rotation. The findings support the development of DE imaging capability for CBCT of the extremities in areas conventionally in the domain of MRI.

High-fidelity artifact correction for cone-beam CT imaging of the brain.

CT is the frontline imaging modality for diagnosis of acute traumatic brain injury (TBI), involving the detection of fresh blood in the brain (contrast of 30-50 HU, detail size down to 1 mm) in a non-contrast-enhanced exam. A dedicated point-of-care imaging system based on cone-beam CT (CBCT) could benefit early detection of TBI and improve direction to appropriate therapy. However, flat-panel detector (FPD) CBCT is challenged by artifacts that degrade contrast resolution and limit application in soft-tissue imaging. We present and evaluate a fairly comprehensive framework for artifact correction to enable soft-tissue brain imaging with FPD CBCT. The framework includes a fast Monte Carlo (MC)-based scatter estimation method complemented by corrections for detector lag, veiling glare, and beam hardening.The fast MC scatter estimation combines GPU acceleration, variance reduction, and simulation with a low number of photon histories and reduced number of projection angles (sparse MC) augmented by kernel de-noising to yield a runtime of ~4 min per scan. Scatter correction is combined with two-pass beam hardening correction. Detector lag correction is based on temporal deconvolution of the measured lag response function. The effects of detector veiling glare are reduced by deconvolution of the glare response function representing the long range tails of the detector point-spread function. The performance of the correction framework is quantified in experiments using a realistic head phantom on a testbench for FPD CBCT.Uncorrected reconstructions were non-diagnostic for soft-tissue imaging tasks in the brain. After processing with the artifact correction framework, image uniformity was substantially improved, and artifacts were reduced to a level that enabled visualization of ~3 mm simulated bleeds throughout the brain. Non-uniformity (cupping) was reduced by a factor of 5, and contrast of simulated bleeds was improved from ~7 to 49.7 HU, in good agreement with the nominal blood contrast of 50 HU. Although noise was amplified by the corrections, the contrast-to-noise ratio (CNR) of simulated bleeds was improved by nearly a factor of 3.5 (CNR = 0.54 without corrections and 1.91 after correction). The resulting image quality motivates further development and translation of the FPD-CBCT system for imaging of acute TBI.

Assessment of image quality in soft tissue and bone visualization tasks for a dedicated extremity cone-beam CT system.

OBJECTIVE: To assess visualization tasks using cone-beam CT (CBCT) compared to multi-detector CT (MDCT) for musculoskeletal extremity imaging. METHODS: Ten cadaveric hands and ten knees were examined using a dedicated CBCT prototype and a clinical multi-detector CT using nominal protocols (80 kVp-108mAs for CBCT; 120 kVp- 300 mAs for MDCT). Soft tissue and bone visualization tasks were assessed by four radiologists using five-point satisfaction (for CBCT and MDCT individually) and five-point preference (side-by-side CBCT versus MDCT image quality comparison) rating tests. Ratings were analyzed using Kruskal-Wallis and Wilcoxon signed-rank tests, and observer agreement was assessed using the Kappa-statistic. RESULTS: Knee CBCT images were rated "excellent" or "good" (median scores 5 and 4) for "bone" and "soft tissue" visualization tasks. Hand CBCT images were rated "excellent" or "adequate" (median scores 5 and 3) for "bone" and "soft tissue" visualization tasks. Preference tests rated CBCT equivalent or superior to MDCT for bone visualization and favoured the MDCT for soft tissue visualization tasks. Intraobserver agreement for CBCT satisfaction tests was fair to almost perfect (κ ~ 0.26-0.92), and interobserver agreement was fair to moderate (κ ~ 0.27-0.54). CONCLUSION: CBCT provided excellent image quality for bone visualization and adequate image quality for soft tissue visualization tasks. KEY POINTS: • CBCT provided adequate image quality for diagnostic tasks in extremity imaging. • CBCT images were "excellent" for "bone" and "good/adequate" for "soft tissue" visualization tasks. • CBCT image quality was equivalent/superior to MDCT for bone visualization tasks.

High-Performance Soft-Tissue Imaging in Extremity Cone-Beam CT.

PURPOSE: Clinical performance studies of an extremity cone-beam CT (CBCT) system indicate excellent bone visualization, but point to the need for improvement of soft-tissue image quality. To this end, a rapid Monte Carlo (MC) scatter correction is proposed, and Penalized Likelihood (PL) reconstruction is evaluated for noise management. METHODS: The accelerated MC scatter correction involved fast MC simulation with low number of photons implemented on a GPU (107 photons/sec), followed by Gaussian kernel smoothing in the detector plane and across projection angles. PL reconstructions were investigated for reduction of imaging dose for projections acquired at ~2 mGy. RESULTS: The rapid scatter estimation yielded root-mean-squared-errors of scatter projections of ~15% of peak scatter intensity for 5·106 photons/projection (runtime ~0.5 sec/projection) and 25% improvement in fat-muscle contrast in reconstructions of a cadaveric knee. PL reconstruction largely restored soft-tissue visualization at 2 mGy dose to that of 10 mGy FBP image. CONCLUSION: The combination of rapid (5-10 minutes/scan) MC-based, patient-specific scatter correction and PL reconstruction offers an important means to overcome the current limitations of extremity CBCT in soft-tissue imaging.

Peripheral Quantitative CT (pQCT) Using a Dedicated Extremity Cone-Beam CT Scanner.

PURPOSE: We describe the initial assessment of the peripheral quantitative CT (pQCT) imaging capabilities of a cone-beam CT (CBCT) scanner dedicated to musculoskeletal extremity imaging. The aim is to accurately measure and quantify bone and joint morphology using information automatically acquired with each CBCT scan, thereby reducing the need for a separate pQCT exam. METHODS: A prototype CBCT scanner providing isotropic, sub-millimeter spatial resolution and soft-tissue contrast resolution comparable or superior to standard multi-detector CT (MDCT) has been developed for extremity imaging, including the capability for weight-bearing exams and multi-mode (radiography, fluoroscopy, and volumetric) imaging. Assessment of pQCT performance included measurement of bone mineral density (BMD), morphometric parameters of subchondral bone architecture, and joint space analysis. Measurements employed phantoms, cadavers, and patients from an ongoing pilot study imaged with the CBCT prototype (at various acquisition, calibration, and reconstruction techniques) in comparison to MDCT (using pQCT protocols for analysis of BMD) and micro-CT (for analysis of subchondral morphometry). RESULTS: The CBCT extremity scanner yielded BMD measurement within ±2-3% error in both phantom studies and cadaver extremity specimens. Subchondral bone architecture (bone volume fraction, trabecular thickness, degree of anisotropy, and structure model index) exhibited good correlation with gold standard micro-CT (error ~5%), surpassing the conventional limitations of spatial resolution in clinical MDCT scanners. Joint space analysis demonstrated the potential for sensitive 3D joint space mapping beyond that of qualitative radiographic scores in application to non-weight-bearing versus weight-bearing lower extremities and assessment of phalangeal joint space integrity in the upper extremities. CONCLUSION: The CBCT extremity scanner demonstrated promising initial results in accurate pQCT analysis from images acquired with each CBCT scan. Future studies will include improved x-ray scatter correction and image reconstruction techniques to further improve accuracy and to correlate pQCT metrics with known pathology.

SU-E-I-101: Initial Implementation and Evaluation of AAPM TG-150 Draft Image Receptor Non-Uniformity Testing Recommendations.

PURPOSE: To implement in software the procedures described in AAPM Task Group 150's draft recommendations for image receptor performance testing, and to evaluate the effectiveness and practicality of these procedures. METHODS: Images of flat fields were acquired using digital x-ray image receptors at 6 cooperating institutions. Four flat field images obtained with each detector spanned a range of input detector air kerma. Software based on AAPM TG150's draft report processed the test images and generated results. Image receptor response and several measures of non-uniformity were evaluated. Images were divided into 10 mm square regions, after eliminating 10 mm borders. For each region, signal (mean), noise (standard deviation) and SNR were calculated. Characteristic signal, noise and SNR were calculated based on average values from all regions. Local non-uniformity for signal (SLN), noise (NLN) and SNR (SNRLN) were expressed as the maximum ratio of the absolute difference between each region's value and its 4 nearest neighbors, to the respective characteristic value. Global non-uniformity (SGN, NGN, SNRGN) were expressed similarly but differences between maximum and minimum values obtained from the regions were used (without comparison to local neighbors). RESULTS: TG150 tests discriminated between good and poorly performing detectors. Improper detector calibration was detectable, with noise non-uniformity proving to be a more sensitive measure than signal or SNR non-uniformity. Detector rotation relative to calibration conditions produced a greater change in signal non-uniformity than the other measures. Image receptor structured noise was characterized by an increase in noise non-uniformity with incident air kerma. CONCLUSIONS: AAPM TG150's proposed approach to image receptor testing was implemented and evaluated. The approach appears to be an effective and practical one for routine quality assurance testing of digital radiographic image receptors.

SU-E-I-103: Detecting Anomalous Pixels and Correlated Artifacts in Digital Detectors from Flat-Field Images.

PURPOSE: Anomalous pixels may be defined as those pixels whose exposure response relationship is deviant from the typical, expected or calibrated response. A group of anomalous pixels may Result in visible correlated artifacts. Here we demonstrate an approach to identify anomalous pixels and correlated artifacts using flat-field images. METHODS: Using manufacturer specific calibration geometry, sets of four flat-field images per detector were obtained with varying input air kerma values (0.5 to 160 µGy) from 9 digital detectors at 6 institutions. Images obtained before and after calibration, with both proper and improper gain maps and structured artifacts were additionally acquired with some detectors. Image analysis methodology under consideration by AAPM Task Group 150 was used.After eliminating 10mm borders, images were divided into square regions (100mm2 ). Anomalous pixels were identified as pixels within each region with valuesabove or below ±3 standard deviations (SD) relative to the mean value of the region. If these pixels were identified in all four images comprising a set, then they were reported as anomalous. Line artifacts were identified as rows and columns with cumulative profile values that were above or below ±3 SD with respect to the mean value of neighboring profiles in the set of four flat-field mages. Results were verified with visual inspection of the images. RESULTS: For four sets of images, the algorithm did not identify any anomalous pixels, and none were spotted on visible inspection as well, while for five sets of images the identified anomalous pixels matched visual inspection results. Anomalous pixel detection failed in regions with an unusually large number of defects and structured noise, since those regions exhibited relatively large SD. Line artifacts consistent with visual analysis were identified correctly when present. CONCLUSIONS: A practical approach to identify anomalous pixels and correlated artifacts from flat-field images is demonstrated.

WE-G-217BCD-04: Diagnostic Image Quality Evaluation of a Dedicated Extremity Cone- Beam CT Scanner: Pre-Clinical Studies and First Clinical Results.

PURPOSE: To assess the diagnostic performance of a prototype cone-beam CT (CBCT) scanner developed for musculoskeletal extremity imaging. Studies involved controlled observer studies conducted subsequent to rigorous technical assessment as well as patient images from the first clinical trial in imaging the hand and knee. METHODS: Performance assessment included: 1.) rigorous technical assessment; 2.) controlled observer studies using CBCT images of cadaveric specimens; and 3.) first clinical images. Technical assessment included measurement of spatial resolution (MTF), constrast, and noise (SDNR) versus kVp and dose using standard CT phantoms. Diagnostic performance in comparison to multi- detector CT (MDCT) was assessed in controlled observer studies involving 12 cadaveric hands and knees scanned with and without abnormality (fracture). Observer studies involved five radiologists rating pertinent diagnostics tasks in 9-point preference and 10-point diagnostic satisfaction scales. Finally, the first clinical images from an ongoing pilot study were assessed in terms of diagnostic utility in disease assessment and overall workflow in patient setup. RESULTS: Quantitative assessment demonstrated sub-mm spatial resolution (MTF exceeding 10% out to 15-20 cm-1) and SDNR sufficient for relevant soft-tissue visualization tasks at dose <10 mGy. Observer studies confirmed optimal acquisition techniques and demonstrated superior utility of combined soft-tissue visualization and isotropic spatial resolution in diagnostic tasks. Images from the patient trial demonstrate exquisite contrast and detail and the ability to detect tissue impingement in weight-bearing exams. CONCLUSIONS: The prototype CBCT scanner provides isotropic spatial resolution superior to standard-protocol MDCT with soft-tissue visibility sufficient for a broad range of diagnostic tasks in musculoskeletal radiology. Dosimetry and workflow were advantageous in comparison to whole-body MDCT. Multi-mode and weight-bearing capabilities add valuable functionality. An ongoing clinical study further assesses diagnostic utility and defines the role of such technology in the diagnostic arsenal. - Research Grant, Carestream Health - Research Grant, National Institutes of Health 2R01-CA-112163.

Task-based modeling and optimization of a cone-beam CT scanner for musculoskeletal imaging.

PURPOSE: This work applies a cascaded systems model for cone-beam CT imaging performance to the design and optimization of a system for musculoskeletal extremity imaging. The model provides a quantitative guide to the selection of system geometry, source and detector components, acquisition techniques, and reconstruction parameters. METHODS: The model is based on cascaded systems analysis of the 3D noise-power spectrum (NPS) and noise-equivalent quanta (NEQ) combined with factors of system geometry (magnification, focal spot size, and scatter-to-primary ratio) and anatomical background clutter. The model was extended to task-based analysis of detectability index (d') for tasks ranging in contrast and frequency content, and d' was computed as a function of system magnification, detector pixel size, focal spot size, kVp, dose, electronic noise, voxel size, and reconstruction filter to examine trade-offs and optima among such factors in multivariate analysis. The model was tested quantitatively versus the measured NPS and qualitatively in cadaver images as a function of kVp, dose, pixel size, and reconstruction filter under conditions corresponding to the proposed scanner. RESULTS: The analysis quantified trade-offs among factors of spatial resolution, noise, and dose. System magnification (M) was a critical design parameter with strong effect on spatial resolution, dose, and x-ray scatter, and a fairly robust optimum was identified at M ∼ 1.3 for the imaging tasks considered. The results suggested kVp selection in the range of ∼65-90 kVp, the lower end (65 kVp) maximizing subject contrast and the upper end maximizing NEQ (90 kVp). The analysis quantified fairly intuitive results-e.g., ∼0.1-0.2 mm pixel size (and a sharp reconstruction filter) optimal for high-frequency tasks (bone detail) compared to ∼0.4 mm pixel size (and a smooth reconstruction filter) for low-frequency (soft-tissue) tasks. This result suggests a specific protocol for 1 × 1 (full-resolution) projection data acquisition followed by full-resolution reconstruction with a sharp filter for high-frequency tasks along with 2 × 2 binning reconstruction with a smooth filter for low-frequency tasks. The analysis guided selection of specific source and detector components implemented on the proposed scanner. The analysis also quantified the potential benefits and points of diminishing return in focal spot size, reduced electronic noise, finer detector pixels, and low-dose limits of detectability. Theoretical results agreed quantitatively with the measured NPS and qualitatively with evaluation of cadaver images by a musculoskeletal radiologist. CONCLUSIONS: A fairly comprehensive model for 3D imaging performance in cone-beam CT combines factors of quantum noise, system geometry, anatomical background, and imaging task. The analysis provided a valuable, quantitative guide to design, optimization, and technique selection for a musculoskeletal extremities imaging system under development.

A dedicated cone-beam CT system for musculoskeletal extremities imaging: design, optimization, and initial performance characterization.

PURPOSE: This paper reports on the design and initial imaging performance of a dedicated cone-beam CT (CBCT) system for musculoskeletal (MSK) extremities. The system complements conventional CT and MR and offers a variety of potential clinical and logistical advantages that are likely to be of benefit to diagnosis, treatment planning, and assessment of therapy response in MSK radiology, orthopaedic surgery, and rheumatology. METHODS: The scanner design incorporated a host of clinical requirements (e.g., ability to scan the weight-bearing knee in a natural stance) and was guided by theoretical and experimental analysis of image quality and dose. Such criteria identified the following basic scanner components and system configuration: a flat-panel detector (FPD, Varian 3030+, 0.194 mm pixels); and a low-power, fixed anode x-ray source with 0.5 mm focal spot (SourceRay XRS-125-7K-P, 0.875 kW) mounted on a retractable C-arm allowing for two scanning orientations with the capability for side entry, viz. a standing configuration for imaging of weight-bearing lower extremities and a sitting configuration for imaging of tensioned upper extremity and unloaded lower extremity. Theoretical modeling employed cascaded systems analysis of modulation transfer function (MTF) and detective quantum efficiency (DQE) computed as a function of system geometry, kVp and filtration, dose, source power, etc. Physical experimentation utilized an imaging bench simulating the scanner geometry for verification of theoretical results and investigation of other factors, such as antiscatter grid selection and 3D image quality in phantom and cadaver, including qualitative comparison to conventional CT. RESULTS: Theoretical modeling and benchtop experimentation confirmed the basic suitability of the FPD and x-ray source mentioned above. Clinical requirements combined with analysis of MTF and DQE yielded the following system geometry: a -55 cm source-to-detector distance; 1.3 magnification; a 20 cm diameter bore (20 x 20 x 20 cm3 field of view); total acquisition arc of -240 degrees. The system MTF declines to 50% at -1.3 mm(-1) and to 10% at -2.7 mm(-1), consistent with sub-millimeter spatial resolution. Analysis of DQE suggested a nominal technique of 90 kVp (+0.3 mm Cu added filtration) to provide high imaging performance from -500 projections at less than -0.5 kW power, implying -6.4 mGy (0.064 mSv) for low-dose protocols and -15 mGy (0.15 mSv) for high-quality protocols. The experimental studies show improved image uniformity and contrast-to-noise ratio (without increase in dose) through incorporation of a custom 10:1 GR antiscatter grid. Cadaver images demonstrate exquisite bone detail, visualization of articular morphology, and soft-tissue visibility comparable to diagnostic CT (10-20 HU contrast resolution). CONCLUSIONS: The results indicate that the proposed system will deliver volumetric images of the extremities with soft-tissue contrast resolution comparable to diagnostic CT and improved spatial resolution at potentially reduced dose. Cascaded systems analysis provided a useful basis for system design and optimization without costly repeated experimentation. A combined process of design specification, image quality analysis, clinical feedback, and revision yielded a prototype that is now awaiting clinical pilot studies. Potential advantages of the proposed system include reduced space and cost, imaging of load-bearing extremities, and combined volumetric imaging with real-time fluoroscopy and digital radiography.

Multiscale deformable registration for dual-energy x-ray imaging.

Dual-energy (DE) imaging of the chest improves the conspicuity of subtle lung nodules through the removal of overlying anatomical noise. Recent work has shown double-shot DE imaging (i.e., successive acquisition of low- and high-energy projections) to provide detective quantum efficiency, spectral separation (and therefore contrast), and radiation dose superior to single-shot DE imaging configurations (e.g., with a CR cassette). However, the temporal separation between high-energy (HE) and low-energy (LE) image acquisition can result in motion artifacts in the DE images, reducing image quality and diminishing diagnostic performance. This has motivated the development of a deformable registration technique that aligns the HE image onto the LE image before DE decomposition. The algorithm reported here operates in multiple passes at progressively smaller scales and increasing resolution. The first pass addresses large-scale motion by means of mutual information optimization, while successive passes (2-4) correct misregistration at finer scales by means of normalized cross correlation. Evaluation of registration performance in 129 patients imaged using an experimental DE imaging prototype demonstrated a statistically significant improvement in image alignment. Specific to the cardiac region, the registration algorithm was found to outperform a simple cardiac-gating system designed to trigger both HE and LE exposures during diastole. Modulation transfer function (MTF) analysis reveals additional advantages in DE image quality in terms of noise reduction and edge enhancement. This algorithm could offer an important tool in enhancing DE image quality and potentially improving diagnostic performance.

Cardiac gating with a pulse oximeter for dual-energy imaging.

The development and evaluation of a prototype cardiac gating system for double-shot dual-energy (DE) imaging is described. By acquiring both low- and high-kVp images during the resting phase of the cardiac cycle (diastole), heart misalignment between images can be reduced, thereby decreasing the magnitude of cardiac motion artifacts. For this initial implementation, a fingertip pulse oximeter was employed to measure the peripheral pulse waveform ('plethysmogram'), offering potential logistic, cost and workflow advantages compared to an electrocardiogram. A gating method was developed that accommodates temporal delays due to physiological pulse propagation, oximeter waveform processing and the imaging system (software, filter-wheel, anti-scatter Bucky-grid and flat-panel detector). Modeling the diastolic period allowed the calculation of an implemented delay, t(imp), required to trigger correctly during diastole at any patient heart rate (HR). The model suggests a triggering scheme characterized by two HR regimes, separated by a threshold, HR(thresh). For rates at or below HR(thresh), sufficient time exists to expose on the same heartbeat as the plethysmogram pulse [t(imp)(HR) = 0]. Above HR(thresh), a characteristic t(imp)(HR) delays exposure to the subsequent heartbeat, accounting for all fixed and variable system delays. Performance was evaluated in terms of accuracy and precision of diastole-trigger coincidence and quantitative evaluation of artifact severity in gated and ungated DE images. Initial implementation indicated 85% accuracy in diastole-trigger coincidence. Through the identification of an improved HR estimation method (modified temporal smoothing of the oximeter waveform), trigger accuracy of 100% could be achieved with improved precision. To quantify the effect of the gating system on DE image quality, human observer tests were conducted to measure the magnitude of cardiac artifact under conditions of successful and unsuccessful diastolic gating. Six observers independently measured the artifact in 111 patient DE images. The data indicate that successful diastolic gating results in a statistically significant reduction (p < 0.001) in the magnitude of cardiac motion artifact, with residual artifact attributed primarily to gross patient motion.